function TFN=calcTFNAB(d1,d2,d3,d4)
T=[1 2 3 4;2 5 6 7;3 6 8 9;4 7 9 10];
A=T(d1,d3);
B=T(d2,d4);
K=[  1     2     3     4     5     6     7     8     9    10    11    11    13    12    13    14    17    15    19  16    17    22    23    18    19    26    20    21    29    22    31    23    33    34    24    25    37    38   39    26    41    27    43    44    28    46    47    29    30    31    51    52    53    32    55    33    57   58    59    34    61    62    35    36    65    66    67    68    69    37    71    38    73    74    75    76    77    78    79    39    40    82    83    84    85    86    87    88    89    41    91    92    93    94    95    96    97    98    99    42];
TFN=K(A*B);


